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Abstract 

We develop a new approach to construct the operator on lattice for the calcula- 
tion of glueball mass, which is based on the connection between the continuum 
limit of the chosen operator and the quantum number J^^ of the state stud- 
ied. The spin of the state studied in this approach is then determined uniquely 
and directly in numerical simulation. Furthermore, the approach can be ap- 
plied to calculate the mass of glueball states (ground or excited states) with 
any spin J including J > 4. Under the quenched approximation, we present 
a pre-calculation result for the mass of 0"'""'" state and 2"'""'" state, which are 
1754(85) (86)Me\/ and 2417(56) (117)MeV^, respectively. 



1 Introduction 



During the past two decades, extensive Monte Carlo simulations were carried out to 
calculate the glueball mass spectra on lattice. B B B Most of these papers car- 
ried out two key steps: one is the choice of glueball operators with certain quantum 
number J^*^ basing on the method introduced in Ref. and the other is the appli- 
cation of variational principle. Meanwhile, with a great amount of the improvement, 
such as fuzzying and smearing, etc., these approaches surely work well and many 
results of well-controlled errors are obtained. However, the simulations are still puz- 
zled by an ambiguity, i.e., how to identify the definite spin J from the irreducible 
representation R of the cubic point group, according to which the glueball operators 
transform, because the corresponding between this irreducible representation R of 
cubic point group and spin J is not one-to-one. 

Meanwhile, basing on the representation theory of 0(4) group and the hypercu- 
bic group, Mandula et. al.0 developed an elegant scheme for the choice of glueball 
operators. With the definition of lattice gauge field, they utilized the lattice color 
electric and magnetic field to construct operators with definite J^*" through the 
decomposition of the composition of lattice color electric and magnetic field into 
certain representation of hypercubic group. But, the correspondence between irre- 
ducible representation of hypercubic group and spin J is also not one-to-one. The 
'leading spin' is then assumed when a — 0. But, this assumption cannot deter- 
mine the spin uniquely. For example, one don't know how to separate 'leading spin' 
J = 1~ from 'leading spin' J = 2"*" in 6*^"'"'' representation and one can not get the 
content of the non-leading spinO. 

We would like to show a possible solution to these troubles in this paper. Unlike 
above references, we start our discussion with the asymptotic expansion of the op- 
erator. By expanding the chosen operator according to power of lattice spacing a, 
we require that the leading term of the expansion of the chosen operator belongs to 
the irreducible representation J^^ of S'0(3)^*^ group. We assume that the leading 
term of the expansion will give the main effect to the state studied when a — 0, 
and the contribution should be only given by the leading term in the continuum. 
Therefore, the spin of the corresponding state is uniquely determined by the leading 
term of the expansion when the lattice tends to continuum. 

We claim that there are two advantages in this approach: (1) Appending to the 
application of variational method one can also study the contribution of the different 
current ( with the same J^^) to the definite state; (2) One can determine the spin of 
the corresponding state unambiguously and directly. In this way, we can distinguish 
state with definite J^^ from other J^^ states. 

Some observations are shown in the forthcoming section. We introduce our 
method in section 3 and give an example of pre-calculation to verify our method in 
section 4. Section 5 is a short summary. 



2 Some Observations 



In the continuum, the glueball states with definite quantum number J^^ make up 
of the bases of certain irreducible representation J^'^ of SO{dt)^^ group. But, on 
lattice, there is only its finite point subgroup, O^*", and its corresponding irre- 
ducible representations R^'^{R = Ai, A2, E, Ti and T2). Then to measure glueball 
mass in lattice QCD, there arises such problem: how we get correct results by 
only utilizing O^'" group. To solve the problem, authors make the continuum limit 
assumption(/5 = 00) since Berg and BilloireS, i.e., denoting masses of the states 
extracted from operators in the irreducible representation R by m{R) and masses 
of the states with certain spin J by m{J), they assume: |^, |^ 

m(0^^) = m(Af^), (la) 

m(l^^) = m{T[^), (16) 

m(2^^) = m(E^^) = m(T/'^), (Ic) 

m(3^^) = m(Af^). (Id) 

But as Ref. shows, this assumption is not always right. From simulation results, 
in Ref. 0, for example, T^^ channel was not interpreted as J^'" = 1+"^ states but 
most likely as J^'-' = 3"*"^ state ( less likely J = 6, 7, 9, ... interpretation cannot rule 
out ), since it seems this channel degenerated with A2"^ channel in the continuum. 
As to the Ti~~ channel, it was interpreted as J^*" = l"*"" state instead of 3"*"". But, 
A2~,Ti * and channels become degenereted in the continuum and they were 
interpreted as 3"''" state. Therefore, we hope to develop an approach to determine 
spin J^*" of states studied directly and uniquely in lattice simulation. 

Meanwhile, on a D = 2 + 1 lattice, Johnson and Teper0 found that in Af^ 
channel there exist two states with different masses. They interpret the higher one 
as 4+"*" state and lower one as 0+"^ state. Therefore, they also consider that one 
needs a systematic and general procedure to construct operators of arbitrary spin 
as a ^ olll. 



Now, we present a possible procedure to solve these problems here. Let us begin 
the discussion with some observations. 



A. 

An arbitrary state {ip > can be regarded as generated by current o acting on vacuum 
|0 >: 

1^ >= o|0 > . (2) 

Since |0 > is invariance under Poincare group and SUc{S) group, the character 
of IV' > can be described by o. For simplicity, we only consider currents with 
mass dimension 4 here, saying Bf{x)Bj{x), where Bf is color magnetic field and 
i, j = 1, 2, 3. 



Since \ip > is color singlet, we also require o is color singlet. One can get 6 color 
singlet currents from above combinations: 

2Tr{B,B,) = Y: BtB], 

a=l 

where Bi = E 5f ^. 

a=l 

We also require \ip > and o transform as certain representation J^*^ under 
SO (3)^'-' group. Since B transforms as 1"*"" under this group, The P, C of Tr{BiBj) 
are ++. Using C-G coefficients, we can decompose these bases consisting of Tr{BiBj) 
into J = and J = 2( Due to color singlet, there is no bases with J = 1). Q Then, 
in the subduced representation J I O of the rotation group 5*0(3) restricted to sub- 
group O, we find that the base in J = is the basis of representation Ai, and we can 
further reduce another five bases in J = 2 according to irreducible representations 
E and T2 of the cubic point group. So, we can categorize these bases as: 

J = 0: au = Tr{BiBi + B2B2 + BsB;)] (3a) 
J = 2: ei = Tr{B,Bi - B2B2), 62 = Tr{BiBi + B2B2 - 2B^B^)- (36) 
t2i = Tr{B2Bs), t22 = TriBiBs), taa = Tr{B^B2). (3c) 

Here an is the basis of representation Ai, e\ and 62 construct bases of representation 
E, while ^21, ^22 and ^23 just make up of bases of representation T2. 

This is just what table 3 ( see Appendix) tells us: the subduced representation 
J = 2 of the rotation group can be decomposing into representation E and T2 in 
the cubic group; the subduced representation J = is just representation Ai. 

We can make similar analysis for higher mass-dimensional gauge invariant oper- 
ator consisting of color magnetic field and its covariant derivative. 

B. 

Now, we consider how to construct a glueball operator on lattice. By expanding the 
chosen operator according to power of spacing a, we require that the leading term of 
the expansion of the chosen operator belongs to and only belongs to the irreducible 
representation J^^ of 5*0(3)^'" group. This is the key point of this paper. Therefore, 
we first consider the perturbative expansion of Wilson loops on lattice according to 
power of lattice spacing a. A simple example is plaquette operator (We denote unit 
vector in the positive i-direction by i): 

O^, = EO.,H = E^^[l - V{n,i)\J{n +?,j)f/-i(n + j,z)f/-i(n,j)], (4) 

n n 

where the link variable f/ is a connector and defined by 

/■a 

U{n,i) = P exp{i / dtAi{an + it), (5) 
Jo 

^ The systemic decomposition was well discussed by Jaffe et al^^ in the study of the qualitative 
features of the glueball spectrum. They suggest to construct glueball operators for certain J^'-' 
states with color magnetic and electric fields in the continuum case. But, we should go further to 
study the construction of operators on lattice as the rest of the paper points out. 



Where P is path-order operator. 

There are two methods to expand the operator. One is the apphcation of non- 
Abehan Stokes theorem^ S and the other one is introduced by Luscher and Weisz 
in Ref. 0. We use both methods and get the same resuhs: 

0^, = E{^Tr(F,,F,,)H + f Tr(F,,F,,F,,)(n) 

n 

-iTr{F,M + D])F,,){n)} + 0(a«), (6) 

where Fij = diAj — djAi — i[Ai, Aj] is gauge field and Di- = di — i[Ai, ■] is covariant 
derivative. 

We now consider the PC = ++ sector of operators, or real part of operators 
in Eq. (^) with ignoring the second term of r.h.s. in Eq. (|^). Due to Oji = 0*j, 
there are three non-zero independent operators ReOu, Re023, ReO^i. Restricting 
oneself into the cubic group, one can combine these operators into representation 
Af^ and 

A++ : Re{023 + Oi3 + O23) = ^ETr{B,Bi + B2B2 + ^gfis)!^) + 0{a^y, 

n 

: Rei023 - O13) = ^ ETr(5i5i - B2B2){n) + 0{a^), 

n 

Re{023 + On - 2O12) = ^ETr{BiB, + ^2^2 - 2B3B3){n) + 0{a'), (7) 

n 

where color magnetic field is Bi = ^ijkPjk- 

jk 

We suppose that the leading term gives the most contribution of the operator 
when a is small enough, or, only the leading term gives the contribution in the 
continuum. While comparing Eq. (0) to Eq. (3), it is assured that, in the continuum 
limit, the state extracted from such operator Af~^ is J = 0, and the states extracted 
from the operator in E'^^ corresponds to J = 2 state. 

We should emphasis again, in the general case, the continuum limits of the 
operator in representation E or in T2 is not always corresponding to J = 2, i.e., 
the parallelism in Eq. (1) does not always hold. Only after expanding the chosen 
operator as we do above, we are then able to affirm or disaffirm the parallelism. 

By the way, we should point out here that the non-leading terms in the expansion 
of the operator do not always belong to the same J^*^ as that of leading term, which 
will bring up the mixing with different spin J. But, as argued above, we expect that 
this artificial mixing will decrease with the decreasing of lattice spacing a so that 
the mixing should vanish when a — although it will affect our error estimate. On 
the other hand, we can utilize non-leading terms to explore high-spin states. 

These two examples tell us that to calculate the mass of the definite J^*" state, 
we should require the continuum limits of our operators belong to and only belong 
to J^'" representation of S0{3)^'" group. One can get this aim by using the com- 
bination of the different operators which belong to the same R^'^ . We will present 
an example to construct the operator in the following section and then show the 
simulation results and discuss the errors in section 4. 



3 The Construction of the Operator 

We exemplify here how to construct operator 0^^ and 2"'"^ up to a^. On lattice, 
gauge-invariant current o with 0++ corresponding to the scalar glueball can be writ- 
ten as 

3 

o — '^{a^'^Tr{BiBj){n) -\- x {current with mass dimension 6) -|- • • •}, (8) 

n i=l 

3 

where the current with mass dimension 6 is the combination of J2Tr(DiFijDiFij), 

3 3 

Tr{DiFjkDiFjk) and J2 Tr{DiFikDjFjk). For simphcity, we only consider the 

current o up to mass dimension 4 in this paper. Then, let us observe the sum of the 
planar special 2x1 rectangular over all lattice sites: 

n 

= ^Y.Tr{[l-U{n,i)U{n + i,i)U{n + 2i,j)U-\n + i+j,i) 

U~\n + J, i)U-\n, j)] + [l-U (n, i)U{n + z, j)U{n + i + j, j) 
U-\n + 2j,i)U-\n + j,j)U-\n,m. (9) 



The real part of the expansion for the operator is 
R^O[, = j:{^^Tr{F,jF,j){n) 



n 
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Then, we define 



WTr{F,j{D^ + D'AF,j){n)} + 0{a'). (10) 



e,,(n) = Re{0,,{n) - ^0;.(n)). (11) 
Continuum limit of operator Q^j is 

^iJ-i:-rwTr{F.jFij)+0{a% (12) 

n ^ 

Decomposing Qij into Ai~^ according to traditional method, we get the basis of 
representation Ai~^: 

F = 012 + 613 + 623 

= T.^Tr{B,B, + B2B2 + B,B3) + 0{a'). (13) 



Apparently, the quantum number of continuum limit of F is 0"*"^. In other words, 
F transforms as O^"*" under S'0(3)^*" group up to a^. We expect that the symmetry 
of S'0(3) has been restored when a — >• 0, the extracted state should be mainly given 
by the leading term of F, so the extracted state is 0"'""'" one in the continuum limit. 
Operator F is our aimed operator for O"*"^ state. 

We may also choose the bases Gi and G2 of representation E^^ to measure 
tensor glueball mass as follows. The operators and their expansion are: 

Gi = i?e(e23 - 613) = E ^-^Tr{B,B, - B^B^) + 0(a«), (14) 

and 

G2 = i?e(023 + ei3 - 2012) = E -r7:Tr{B,B, + ^2^2 - 2B,B^) + 0{a^). (15) 
According to (3b), they belong to bases of representation 2++ up to 0{a^). 



4 Simulation Results 



Under the quenched approximation, we perform our calculation on anisotropic lattice 
with improved gluonic action as chosen in Ref. 
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where f3 = g is the QCD couple constant, Ug and Ut are mean link renormal- 

ization parameters(we set Ut = I), ^ = as /at is the aspect ratio, and VLsp includes the 
sum over all spatial plaquettes on the lattice, VLtp indicates the temporal plaquettes, 
Vtsr denotes the planar 2x1 spatial rectangular loops and VLgtr refers to the short 
temporal rectangles( one temporal and two spatial links). More detail is given in 
Ref. 0. In each [3 calculation, we set 2800 sweeps to make configurations reach to 
equilibrium and make measurement once after every four sweeps. Our calculation 
spends 80 bins in which there are 70 measurements after reaching equilibrium. The 
method to set the scale used here is introduced by ref \T1 
measured and found they coincide with those in Ref. 0] 
we adopt the data in Ref. [§] as our energy scale. Table 1 shows the simulation 
parameters 



. For each /3, we have 
in the error bound. So 
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Table 1 The glueball simulation parametersl^ 



. Here we assume vq = A10(20)MeV . 



As argued above, we choose operator F to calculate scalar glueball mass and 
operator Gi and G2 to calculate tensor glueball mass. As usual, we calculate 
the average over sample vacua of a correlation C{t) =< 0|o'^(t)o^(0)|0 >, where 



it) 



o{t)— < 0\o{t)\0 > is the vacuum-subtracted form of the chosen operator, 
to determine masses of the corresponding glueball states with the improvements such 
as fuzzying and smearing. At the same time, following the mean field theoryli^l, 
we also replace link variant U by U/ug in the chosen operators due to the tadpole 
correction. After such programmes, we get our results in each j3 which are shown in 
table 2: 
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Table. 2 Glueball energy meat for each j3. 
The numerals in the brackets are the error estimates. 



Now we comment a little on the error estimate. 

First, our action breaks the rotation symmetry to 0(af , a 



i.e., the upper limit 



of the precision in the calculation is 0(a^, a^). Since as argued by many authors, the 
contribution of 0(af) can be ignored, the upper limit of the precision here is o(a^). 

Second, we ignore terms(currents) with mass dimension 6 in Eq. (H). Due to 
dimensional analysis, the contribution of the terms to error should have a square 
mass suppression'Hl, which will make two effects on our mass measurement. One 



is that we should include it in systematic error in the continuum limit, which needs 
further calculation to get its accurate value. Here we simply expect that it is about 
{-^qcd/MY, where we set Aqcd — 250MeV and M is measured mass. The second 
is that it will takes 0{a1) error when 7^ 0. Since it is not statistical error, its 
contribution to error can be fitted by €20^ + c^a^ + ■ ■ ■. 

From the argument and calculated data, we use the formula m{0^~^, a^) = 1.754— 
1.514(a,/ro) + 1.773(a,/ro)2 and m(2++, a,) = 2.417 + 0.783(a,/ro)^ - 0.787(a,/ro)^ 
( unit: GeV) to fit our data. We present our data and fitting curves in Fig. 1. 

The statistical error is 0.076Ge^ for scalar glueball and 0.044^61^ for tensor 
glueball. According to Ref. systematic error is 1 percent ( from aspect ratio 
). But since our method also gives about 2 and 1 percent systematic error for O"*""*" 
and 2+"*" states respectively, the total systematic error is about 2.2 percent( 39Me\^ 
) and 1.4 percent( SAMeV ) respectively. Therefore, the mass of scalar glueball is 
1.754(85)Gey and the mass of tensor glueball is 2.417(56)GeV^. Including the uncer- 
tainty in r^^ = 410(20)Met;, Our final results are: Mg(0++) = 1754(85) (86)Me1/ 
and mG(2++) = 2417(56)(117)Me1/. 



5 Conclusion 



Basing on the connection of the asymptotic expansion of the operator and the quan- 
tum number J^*" of the extracted state, we present a new approach to construct 



operator on lattice for the calculation of the glueball mass, which may solve the 
ambiguity in the simulation. This approach points out that, in general, to calculate 
the mass of definite J^^ glueball states, first one should write out these currents 
which transform as the representation J^^ under the 5*0 (3)^*^ group in continuum 
and decompose them into irreducible representations R^^ of the group O^^ in the 
subduced representation which is obtained by trivially embedding the O^*-^ group 
into the 5*0(3)^'" group; then one should construct corresponding operators which 
belongs to the representation R^*-^ of the O"^*" group on lattice, its continuum limits 
should be those currents mentioned above. 

To verify our approach, we calculate the scalar and tensor glueball mass under the 
quenched approximation in this approach. Since the continuum limit of operator F is 
actual O'^"'", we affirm the mass extracted from the operator F is scalar glueball mass 
and its value is 1754(85) (86)Mey. With the same reason, the mass extracted from 
operator Giii = 1,2) is the tensor glueball one and its value is 2417(56)(117)Mey. 
These results are consistent with those obtained in references 0, ||, H, |T5[. 

Apparently, there is no radical obstacle to prevent us to calculate the mass of 
states with any spin J including J > 4 in this approach. We will make such study 
systematically. We have first calculated the mass of the ground 4"*"+ states in the 
E~^~^ channel with 2+"*" stete under the quenched approximation in this approach. It 
will be shown elsewhere. 

Of course, the operator, its continuum limit transform as J^*" of S0{3)^'^ group, 
is not unique. For example, one can also construct the operator including color 
electric field. With these operators, one can determine their relative weights of 
contribution by variational principle. But, we did not make such treatment here. 

To compare our results with the experiments, we need to calculate it on an 
unquenched lattice. We should also study the mixing between glueball states and 
normal mesons with the same J^*^. Such works should be finished in future. 
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Appendix 

In the cubic group, the subduced representation J of irreducible representation J 
can be decomposed into irreducible representation R. Table 3 shows the multiplicity 
of decomposing J up to J = 6. It is well known that the subduced representation 
with J > 2 are reducible and only up to J = 3 do new irreducible representation of 
the cubic group show up. 
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up to J = 6 for the cubic group. 



Figure Caption 

Figure 1 Masses of scalar and tensor glueball against the lattice spacing (os/ro)^. 
The fitting curve are m(0"'"^,as) = 1.754 — 1.514(as/ro)^ + 1.773(as/ro)^ for scalar 
gluebaU mass and m(2++,a,) = 2.417 + 0.783(a,/ro)2 - 0.787 {as /ro)\unit : GeV) 
for tensor glueball mass. The masses in continuum limit are 1.754(76)Gey and 
2. 417 (4:4:) GeV if we only consider the statistical error. The top data and curve: 
m(4++,a,) = 3.65 - 1.22{as/rof + 2.74:{as/roy{unit : GeV) is for 4++ glueball 
mass. 



